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ABSTRACT 

We present dynamical scaling relations for a homogeneous and representative sample of 
'^SOO massive galaxies, selected only by stellar mass (> 10^° M©) and redshift (0.025 < 
z < 0.05) as part of the ongoing GALEX Arecibo SDSS Survey. We compare baryonic 
TuUy-Fisher (BTF) and Faber-Jackson (BFJ) relations for this sample, and investigate 
how galaxies scatter around the best fits obtained for pruned subsets of disk-dominated 
and bulge-dominated systems. The BFJ relation is significantly less scattered than the 
BTF when the relations are applied to their maximum samples (for the BTF, only 
galaxies with Hi detections), and is not affected by the inclination problems that 
plague the BTF. Disk-dominated, gas-rich galaxies systematically deviate from the 
BFJ relation defined by the spheroids. We demonstrate that by applying a simple 
correction to the stellar velocity dispersions that depends only on the concentration 
index of the galaxy, we are able to bring disks and spheroids onto the same dynamical 
relation — in other words, we obtain a generalized BFJ relation that holds for all 
the galaxies in our sample, regardless of morphology, inclination or gas content, and 
has a scatter smaller than 0.1 dex. We compare the velocity-size relation for the 
three dynamical indicators used in this work, i.e., rotational velocity, observed and 
concentration-corrected stellar dispersion. We find that disks and spheroids are offset 
in the stellar dispersion-size relation, and that the offset is removed when corrected 
dispersions are used instead. The generalized BFJ relation represents a fundamental 
correlation between the global dark matter and baryonic content of galaxies, which is 
obeyed by all (massive) systems regardless of morphology. 

Key words: galaxies: kinematics and dynamics - galaxies: evolution - galaxies: 
fundamental parameters - radio lines: galaxies 



1 INTRODUCTION 

The observed global properties of galaxies obey a diverse set 
of scaling relations, which are fundamental tools to constrain 
models of galaxy formation and evolution. Particularly in- 
* bcatinel@mpa-garching.mpg.de teresting are the correlations between dynamics and lumi- 
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nosity /stellar mass o r size, such as the TuUy-Fisher (TF; 
ITuUv fc Fisheij 1977) relation for s pirals, a nd the Faber- 
Jacks on (FJ; Taber fc Jackson'' 1976^ ), Dn-a dP rcssler et al.' 
1987^■ or fundament al plane (FP; D iorgovski fc Davis 1987: 



Dressier et al] Il987l ) relations for spheroids, because they 



link the luminous to the total mass of galaxies, thus provid- 
ing insights into the interplay between their baryonic and 
dark matter components. 

Historically, the importance of the tight TF and FP 
relations as secondary distance indicators has driven the 
community to assemble galaxy samples that meet strict 
selection criteria, in order to minimize systematic errors 
and scatter. While FP studies targeted elliptical and SO 
galaxies, the samples used for TF applications have typi- 
cally been restricted to late-type spirals with inclinations 
to the line-of-sight larger than 30-40°, preferably observed 
in red or infrared photometric bands to minimize extinc- 
tion effects (e.g.. ICourteaul 1 19971 : iGiovaneUi et all Il997bl : 
[Masters et al. I l2006l ). These data sets are not ideal for 
characterizing the statistical properties of galaxies in gen- 
eral, because the excessive pruning means that they are 
not fair samples of the local Universe. This issue is par- 
ticularly problematic for the comparison with theoretical 
studies and numerical simulations of galaxy formation and 
evolution, which should be based on representative sam- 
ples. For instance, TF studies of different classes of ob- 
jects , such as SOs and early - type spirals (e.g . . iNeistein et al.l 
Il999l : iBedregal et all l2006l : IWiUiams et all I2OI0I. and ref- 
erences therein) , polar ring galaxie s ( lodice et al.l |2003| ) , 



barred spirals ([Courteau et al. 12003), an d gas-rich dwarfs 
l|McGaugh et al.ll2000l : iBegum et al.ll2008l ), have sometimes 



found disagreement with the TF relation of late-type spi- 
rals. Most notably, gas-rich dwarfs lie systematically below 
the TF relation defined by bright galaxies. 

During the last decades, the interest in TF and 
FJ-like scaling relations has shifted from cosmic flow 
applications to constraining galaxy formation models, 
and samples with broader morphological properties have 
been constructed specifically for this purpose (e.g., 
Kannappan et al. ^2002 ; Piza g no et al.l l2007l : iDutton et al.l 
20071 : lAvila-Reese et alJi2008l ). However, current large and 
homogeneous data sets include eithe r spirals (e.g., in addi- 
tion to the works mentioned ab ove, ICourteau et al.l l2007al 
and ISaintonge fc S pckkc"iis 2011*. hereaft er SSll) or early - 
type galaxies (e.g., Bcrnardi et al. , 2003: Graves et al.ir2009l ) 
only. 

There have been attempts to move beyond the spi- 
ral/elliptical dichotomy, and uncover relations between stel- 
lar content and dynam ics that hold for a ll galaxies, inde- 
pendent of morphology. IZaritskv et al. found that all 
galaxies, from disks to spheroids and from dwarf spheroidals 
to giant ellipticals, lie on a two-dimensional surface defined 
by surface brightness, half-light radius, internal velocity and 
mass-to-light ratio. As a measure of internal velocity V, they 
adopt either the rotational velocity Kot for disks or the ve- 
locity dispersion a for spheroids, thus the two types of sys- 
tems are still treated separately (especially because the sam- 
ple is a large but heterogeneous collection of published data 
sets, for which either Kot or a is available). 

It is also important to point out that the search for scal- 
ing relations that are valid for all types of galaxies should 
make use of baryomc masses (i.e., the sum of stellar and 



gas masses) instead of luminosities or stellar masses, be- 
cause they could be more fundamental quantities. This is 
demonstrated by the fact that baryonic scaling relations 
hold for subsets of galaxies that do not follow the corre- 
sponding stellar relations. For example, the offset of the gas- 
rich dwarf galaxies from the stellar TF rel ation disappears 
when their gas mass is t aken into account (iMcGaugh et al.l 
'2000'; Begum et al. 2008'). The baryonic TF relation (BTF; 
iMcGaugh et al. 2000) is linear over 5 orders of magnitude in 
(stellar -|- gas) mass, suggesting that the TF is fundamen- 
tally a relation between baryonic (rather than luminous) 
and total mass of the galaxy. Intriguingly, although sup- 
ported by limited statistics, there is some evidence that gi- 
ant and dwarf ellipticals migh t lie on the same BTF as the 
spirals l|De Riicke et al.ll2007n . Unfortunately, because they 
require estimates of both stellar and gas masses , BTF sam- 
ples ( e.g., .M cGaugh 2005; Geha ct al. 2006; B egum et all 
l2008l : iGurovich et al.l 120041 . l2010i ) are significantly smaUer 
than TF ones. 

To summarize, it is still unclear if disk-dominated galax- 
ies and spheroids obey the same dynamical scaling relations, 
mainly due to the lack of well-defined, representative sam- 
ples of galaxies for which both rotation and stellar disper- 
sion are measured. Certainly, ellipticals might have no gas 
and no detectable rotation, and pure disks might have neg- 
ligible stellar dispersions, but a significant fraction of local 
gala xies (especially mas sive ones) have a disk and a bulge 
(e.g. [Driver et al.l |2007| ). and the dynamical scaling rela- 
tions should account for the smooth transit ion across galajc- 
ies wi th different bulge-to-disk ratios. As ICovington et al.l 
l|2010h point out, a connection between the scaling rela- 
tions of early-type and late-type galaxies is expected on the 
grounds that early-type systems are generally assumed to 
form through mergers of late-type ones. This is especially 
true at the high stellar mass end, where the blue sequence of 
star-forming disks merges onto the red s equence of passively - 
evolving, bulge-dominated galaxies (e.g. iBaldrv et ai1l2006l ). 
and the systems typically host a bulge and a disk. 

In this paper, we investigate dynamical scaling relations 
for a representative sample of ~500 massive galaxies that are 
selected only by stellar mass (M, > 10^'^ M©) and redshift 
(0.025 < z < 0.05), as part of the ongoing GALEX Arecibo 
SDSS Survey (GASS: [Catinella et allboiOl . hereafter Paper 
I). For these galaxies, we have homogeneous measurements 
of structural parameters and ve locity dispersions from the 
Sloan Digital Sky Surv ey (SDSS:lYork et al.l l2000l). NUV-r 
colours firom GALEX (|Martin et al.ll2005l ) and SDSS imag- 
ing, and Hi masses and rotational velocities for the subset 
of objects detected at 21 cm with the Arecibo radio tele- 
scope. This unique sample, which includes massive galaxies 
of all morphological types, allows us to investigate how ob- 
jects that are typically not included in TF or FJ/FP data 
sets scat ter around those relat ions. In the spirit of works lik e 
those of IZaritskv et al.l (|2008l ) and ICovington et all |20lO), 
we wish to establish if there is a fundamental correlation 
between baryonic mass and dynamics that is obeyed by the 
complete galaxy population, regardless of morphology. We 
show that, at least for the massive galaxies in our sample, 
such a relation does exist, and has a scatter smaller than 0.1 
dex, comparable to that of the TF and FJ relations applied 
to their respective pruned subsets. 

This paper is organized as follows. We summarize sam- 
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pie selection and measurements of relevant quantities in § (2] 
We present the baryonic mass-velocity relations, starting 
with the TF and FJ, in § |3l and the velocity-size relations 
in § |4l We discuss our findings and conclude in § [5] 

All the distance-dependent quantities in this work are 
computed assuming Q = 0.3, A = 0.7 and Ho = 70 
km s"^ Mpc"\ 



2 SAMPLE SELECTION AND GALAXY 
PARAMETERS 

The sample used in this work is drawn from GASS, an on- 
going survey which is gathering high-quality Hl-line spectra 
for ~1000 massive galaxies, selected only by stellar mass 
(greater than 10^° Mq) and redshift (0.025 < z < 0.05). 
The GASS targets are located within the intersection of 
the footprints of the SDSS primary spectroscopic survey, 
the projected GALEX Medium Imaging Survey and the 
on-going Hi bhnd Are cibo Legacy Fast ALFA (ALFALFA; 
iGiovanelli et al. 

1 12005!) survey. The galaxies are observed 
with the Arecibo radio telescope until detected or until a 
gcis fraction limit of 1.5 — 5% is reached. For more details 
we refer the reader to Paper I, where the first GASS data 
release (DRl) is presented. 

Here, in addition to the DRl data (176 galaxies), we 
use new Arecibo observations of 240 galaxies that will be 
incorporated in the second GASS data release (Catinella et 
al., in preparation). As discussed in Paper I, GASS does not 
re-observe objects with good Hi detections already avail- 
able from the ALFAL FA survey or the Cornell Hi archive 
l|Springob et al. I l2005l . hereafter S05). To correct the GASS 
sample for its lack of Hl-rich objects, we add galax;ies from 
ALFALFA and S05 in the correct proportions, following the 
procedure detailed in section 7.2 of Paper I. The sample 
thus obtained, which is representative in terms of Hi prop- 
erties, includes 480 galaxies (296 detections and 184 non- 
detections). As explained below, we discard 44 galaxies for 
which the stellar velocity dispersion is not reliable. The fi- 
nal sample includes 436 galaxies (259 detections and 177 
non-detections). 

As mentioned in Paper I, the optical parameters are 
obtained from queries to the SDSS DR7 (|Abazaiian et al] 
I2OO9I ) data base server, unless otherwise noted. Stellar 
masses are derived from SDSS photometry using the spec- 
tra l energy dis tr ibutio n (SED) fittin g technique desc ribed 
in ISalim et al.l (|2007l ). assuming a IChabrierl (|2003l ) ini- 
tial mass function . A v ariety of model SEDs from the 
iBruzual fc CharlotI l|2003 ') library are fitted to each galaxy, 
building a probability distribution for its stellar mass. The 
mean and the width of this distribution are used as mea- 
surements of the stellar mass and its formal error, respec- 
tively. Over the interval probed by GASS, stellar masses are 
believed to be accurate to better than 30 per cent. As a 
measure of galaxy size we adopt -R25, i.e. half the 25 mag 
arcsec"'^ isophote diameter D25 (measured by us on the 
SDSS 3-band images), in kpc. 

The baryonic mass is the sum of stars and gas; the lat- 
ter is computed from the Hi mass, adopting the standard 
1.4 correction factor to account for helium and metals, i.e. 
Afgas ~ 1.4Mhi. This correction neglects the contribution 
of the molecular hydrogen. In our stellar mass regime, the 



amount of H2 does not depend significantly on M* or con- 
centration index and, on ave rage, Mh, /Mhi = 0.3 (but 
with a large scatter, 0.41 dex: ISaintonge et al.ll201ll ). Since 
the contribution of the Hi to the baryonic mass is, for our 
sample, typically small (see §[3TT}, it is safe to neglect the 
H2 (but we did check that including the molecular gas does 
not change our results). We set the Hi masses of the non- 
detections to zero, but we note that using the upper limits 
(which correspond to 1.5-5% of the stellar mass by survey 
design) would make no difference to our plots. Other param- 
eters used in this work are discussed below. 



2.1 Hi line widths 

Hi fine widths from GASS, ALFALFA and the S05 archive 
are measured wi th the same techniqu e, at the 50% of each 
peak level (e.g.. [Catinella et ahlfioOTl , §2.2). However, the 
corrections applied to the raw measurements, as originally 
published, are different. GASS and S05 widths are corrected 
for both instrumental broadening and cosmological redshift, 
whereas ALFALFA widths are correcte d for instrumenta l 
broadening only, following equation 1 in iKent et all (|2008h . 
No turbulent motion or inclination corrections are applied. 
For convenience, we report here the corrections adopted: 

W50 - As 



W, 



GASS,S05 — 



1 + Z 



where W50 is the measured velocity width, z is the galaxy 
redshift, and As is the instrumental broadening correction, 
which differs for the three sources. For GASS Paper I and 
ALFALFA, As is simply the final velocity resolution of the 
spectrum after smoothing (i.e., between 5 and 21 km s~^ 
for GASS spectra, which are Banning and boxcar smoothed, 
and ~10 km s~^ for ALFALFA spectra, which are Hanning 
smoothed only). For S05, As equals 2AiichA, where Avch 
is the channel separation in km s~^ and A is a complex 
function of signal-to-noise ratio (SNR) and type of smooth- 
ing applied. For high SNR, Hanning-smoothed spectra with 
Avch = 5 km s~^ (i.e., the ALFALFA case), A = 0.40, thus 
As is of order of half the velocity resolution after smooth- 
ing. This is smaller than the correction adopted for GASS in 
Paper I, and larger than the ALFALFA one, where the sub- 
traction in quadrature makes the correction negligible (i.e., 
< 1 km s~^) for velocity widths larger than 50 km . 

The issue of homogenizing velocity widths obtained 
from different s ources has been discu ssed in the past, and 
most recently bv lCourtois et aL ||2009| ) for the Extragalactic 
Distance Database ( Tullv et al.l 200 9^) . a large compilation of 
dat a for galaxy distance a nd peculiar motion studies. We fol- 
low [CourtoisIet^L| l|2009l) . who adopt a simplified version of 
the S05 solution, i.e.: As = 2AvX, where A = 0.25 and Aw 
is the final velocity resolution of the spectrum after smooth- 
ing (this is half the correction adopted in Paper I). This is 
in better agreement with our own tests on high SNR GASS 
Hi profiles, where the smoothing was gradually increased 
and the width remeasured. Thus, we keep the S05 widths as 
publishecQ, we adopt the new As correction for the GASS 

^ For the only galaxy in common between SOS and GASS (GASS 
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data, and we uncorrect the ALFALFA widths, H^alfalfaj 
to obtain the raw W50 measurements, then apply the same 
As and {1 + z) corrections adopted for GASS galaxies. 

Lastly, we deproject all velocity widths to edge-on view. 
The inclination to the line-of-sight is computed from: 



150 



(b/a) 



(1) 



where b/a is the minor-to- major axis ratio from the r-band 
exponential fit {expAB.r in the SDSS database), and qo is 
the intrins ic axial ratio of a galaxy seen edge-on. We adop t 
qo = 0.20 (jHolmberelll"958l : see also, e.g. ITuUv et~al]|2009h . 
and set the inclination to 90° for galaxies with b/a < 0.2. 
The value of go = 0.20 applies to galaxies of morphological 
type earlier than Sbc, where as go = 0.13 should be used 
for Sbc and later spirals (e.g. iGiovanelli et al.lll994l . Il997al : 
SSll). However the difference is small, and our choice of go 
seems appropriate for a sample of massive galaxies that was 
not selected for TF studies. 

Thus, rotational velocities are computed as: 



W50 - 0.5Av 
2 (1 + z) sin i ' 



(2) 



where Aw is the final velocity resolution of the spectrum 
after smoothing. 

As noted in Paper I, the Hi masses of galaxies in the 
ALFALFA and SOS archives have been recomputed from the 
tabulated fluxes, in order to be made consistent with GASS 
ones. 

We identified 31 galaxies for which Hi confusion within 
the ~4' Arecibo beam is certain, and therefore the measured 
Hi parameters should not be trusted. Specifically, we closely 
inspected the SDSS images of all the galaxies in our sample, 
and flagged as "confused" those with at least one late-type, 
similar size companion (based on SDSS spectroscopy, i.e. 
with redshift available and within 0.002 of that of the target 
galaxy) within the beam. These objects will be highlighted 
in our analysis when needed. 



2.2 Velocity dispersions 

Galaxy velocity dispersions are measured by fltting stellar 
templates convolved with Gaussian functions to SDSS spec- 
tra, which are obtained through 3"-diameter fiber apertures. 
Because the fiber covers only a fraction of the galaxy light at 
the GASS redshifts, these quantities (catalogued as velDtsp 
in the SDSS spectroscopic data base, and here referred to as 
ffib) need to b e corrected for aperture effects. As comm only 
done (see, e.g. iBernardi et al.|[2003l : [Graves et al.ll2009l . and 
references therein), fiber velocity dispersions are corrected 
to 1/8 effective radius as follows: 



O" = 0"flb 



rah 
ro/S 



(3) 



where r^h ~ 1-5" 

arcseconds, computed as ro = -RdeV vW^O dcv {RdcV and 
(fo/a)dev are the effective radius and axis ratio from the 



and ro is the circular galaxy radius in 

(-RdoV 



38751, AGC 240702), the corrected velocity widtlis agree within 
the quoted errors. 
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Figure 1. Stellar mass distribution for the GASS sample (solid). 
The shaded histogram shows the 44 objects with stellar velocity 
dispersion smaller than 70 km s~^ that are not included in our 
analysis. 



r-band de Vaucouleurs fit, respectively). This correction is 
small, ~3% on average, for the galaxies in our sample. 

Only values of aah above the ~70 km instrumental 
resolution of the SDSS spectra are considered reliable. We 
thus discarded 44 galaxies that do not meet this requirement 
from our sample; the effect of this restriction on the stellar 
mass distribution of the sample is shown in Figure [T] 



3 BARYONIC MASS- VELOCITY RELATIONS 

The main goal of this work is to establish if there is a relation 
between baryonic mass and a measure of the dynamical mass 
(estimated based on rotational velocity, stellar dispersion or 
a combination of the two) that is obeyed by massive galaxies 
regardless of their morphology. 

We begin by comparing the baryonic TF (BTF) and 
baryonic FJ (BFJ) relations for our sample, in order to es- 
tablish which quantity, the Hi rotational velocity or the stel- 
lar dispersion, most reliably traces the baryonic mass of mas- 
sive galaxies. Thanks to its unique selection by stellar mass 
and redshift only, the GASS sample includes massive galax- 
ies of all morphological types, and it is thus ideal for this 
comparison. Indeed, this is not a sample designed for either 
TF or FJ studies, but we can prune it to separately study 
inclined, disk-dominated and spheroidal, bulge-dominated 
systems. Most importantly, we can investigate outliers and 
residuals of the BTF and BFJ relations, which is essential 
in order to determine how to bring disks and spheroids onto 
the same relation. We conclude by illustrating two differ- 
ent ways of obtaining a baryonic relation that holds for all 
the massive galaxies in our sample, and whose dispersion is 
comparable to that of the BTF and BFJ relations. 



3.1 Baryonic Tully-Fisher Relation 

Here we must restrict our sample to galaxies with Hi de- 
tections, for which the rotational velocity can be measured. 
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Figure 2. Relation between baryonic mass and rotational velocity for all the galaxies with Hi detections (a). Magenta diamonds indicate 
objects for which Hi confusion is certain; red and blue stars correspond to GASS 3505 and 35981, respectively (see text). The other 
panels show the same relation for selected subsets: (b) inclined, disk-dominated systems, (c) galaxies with inclinations smaller than 40° 
(cyan) and/or with Vrot < O" (red), and (d) bulge-dominated objects (filled cyan circles are bulge-dominated galaxies with inclination 
smaller than 40°). In the top panels, a dotted (a) or dashed (b) line indicat es the inverse fit to the data points (the scatter is noted in the 
bottom right corner), and the red solid line shows the BTF relation from iMcGaugh et al. I 1I2OOOI ). The dashed line in (b) is reproduced 
in panels (c) and (d) for reference. 



The BTF relation for this sample is shown in Figure [2^. 
We note that, for our interval of stellar masses, the contri- 
bution of the gas to the baryonic mass is generally small 
— if we plotted stellar instead of baryonic masses, most of 
the points would move downward by an amount compara- 
ble to the symbol size (indeed, there are only 13 galaxies in 
our sample with gas fractions Mhi/M*> 50%). We plot the 
BTF as usually done (i.e., baryonic mass versus rotational 
velocity), but we fit the inverse relation, since the scatter is 
clearly associated to the Kot-coordinate (as can be seen by 
comparing Figures[2^ and[5]D, explained below). The inverse 
fit to the data points, i.e. Log Vrot= a Log (A/*+ L4 Mhi) 
+ b, is shown as a dotted linqj. Following SSll, we compute 



^ Tables [Tl and [2l list the expression of the fit, the scatter, and 
the number of galaxies contributing to the sample for all the main 
relations discussed in this work. 



the scatter in the x variable around the best fit by apply- 
ing Tukey's bi-weight, which yields a robust estimate of the 
dispersion in presence of outliers; the scatter (in dex of ve- 
locity) is noted at the bottom of the panel. A lso shown is 
the relation obtained bv lMcGaugh et al] (|2000l ) over nearly 
5 orders of magnitude of baryonic mass (solid line), which 
is in excellent agreement with our data. 

Not surprisingly, the scatter around the best fit is large 
for our data set (0.127 dex), and discarding galaxies with 
confused Hi spectra (magenta symbols; see § 12. ip improves 
it only by a small amount (0.112 dex). However, we recover 
a significantly tighter relation when we prune the sample as 
usually done for TF studies. This is illustrated in Figure 
where a TF subset was obtained by selecting non-confused, 
disk-dominated, inclined galaxies (i.e., objects with concen- 
tration index Rqo/Rso^ 2.8 and inclination i > 40°; the 
scatter of the BTF for this subset is 0.076 dex). We inspected 
the three low- velocity outliers that stand out on the left of 
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the plot. Two of them are galaxies for which the SDSS in- 
clination is clearly overestimated (we measured 17° and 26° 
instead of 46° and 49°, respectively), and the third one has 
uncertain Hi velocity width. Removing these outliers, the 
scatter of the relation becomes 0.072 dex. How does this 
compare with the scatter of other TF samples? There is a 
vast literature on the TF relation, showing that its scatter 
depends on the sample analyzed, velocity indicator, photo- 
metric band, and type of fit performed (forward, inverse, 
bisector or orthogonal). iMcGaugbl (|2005h studied the BTF 
relation using a sample of 60 galaxies with extended Hi ro- 
tation curves, and obtained a scatter of 0.191 dex in solar 
masses from a direct fit (this is the average scatter for the re- 
lations marked as in his Table 2, where the stellar masses 
are computed using stellar population synthesis models as 
we do). For comparison, a direct fit to the data in Figure [2]3 
(excluding the three low- velocity outliers mentioned above) 
yields a scatter of 0.221 dex in solar ma sses. A study that 
is dire ctly comparable to ours is that of I Avila-Reese et al.l 
l|2008l ), who measured a scatter of 0.06 dex in velocity for 
the BTF of a sample of 76 non-interacting disk galaxies, 
with inclinatio ns 35° < i < 80° and fiat rotation curves 
(see also, e.g.. IPizagno et al]|2007l ). Given the crude selec- 
tion of our TF subset (by concentration index), our sample 
is likely to include a broader range of disk galaxy types, thus 
the slightly larger scatter compared to these studies is not 
unexpected. 

Figure [2}; demonstrates that the high- velocity outliers 
of the BTF relation are all galaxies with small inclination 
to the line-of-sight. For these systems, the corrected rota- 
tional velocities become so unreliable that the correlation 
with baryonic mass effectively disappears. Once galaxies 
with inclination smaller than 40° are removed from the sam- 
ple, even bulge-dominated objects lie reasonably close to the 
relation obtained for disk-dominated ones, although with 
larger scatter (empty circles in FigJ2H) . This agrees quali- 
tatively with the findings of IHoI |2003), based on a compi- 
lation of 792 galaxies with inclination larger than 30° and 
heterogeneous measurements from Hyperleda ( iPaturel et al.l 
[2003a , b). Indeed, Ho showed that a ifg-band TF relation ex- 
ists for all Hubble types, including elliptical and SO galaxies 
(although these represent a very small fraction of his sam- 
pie). 

Figure [2}; also shows a population of low- velocity out- 
liers that cannot be explained by small inclinations: these are 
galaxies with rotational velocities that are smaller than the 
stellar velocity dispersions (red symbols). We note that IHoI 
l|2007l ) also identifies a population of galaxies with unusually 
small Vrot/o" ratios, which are outliers in his TF relation. Ex- 
cept for their small rotational velocities compared to their 
central stellar velocity dispersions, his outliers are otherwise 
normal luminous galaxies. These systems might be very in- 
teresting — Ho argues that a significant fraction of their Hi 
gas must be dynamically unrelaxed, having been acquired 
through a minor merger episode or perhaps cold accretion. 
For our sample, the fraction of low Kot/o" outliers is smaller 
(~5%, against 17% for Ho's sample). As shown in the Ap- 
pendix, where these galaxies are described in more detail, 
most of them have asymmetric Hi profiles, suggesting that 
the Hi distribution and/or kinematics might be disturbed 
(although we cannot prove that the gas was externally ac- 
creted) . 



Lastly, here and in other figures of this paper we high- 
light the positions of two interesting galaxies that were dis- 
cussed in our previous works, GASS 3505 and 35981 (marked 
as a red and a blue star in Fig. [2^). These are both unusu- 
ally Hl-rich systems, which are outliers in the gas fraction 
plane discussed in Paper I (relating the Hi mass fraction to 
stellar mass surface density /i* and NUV— r color). However, 
GASS 3505 is an early-type system, whereas GASS 35981 
(UGC 8802) is a disk galaxy with a sharp metallicity drop 
in it s outer disk, which suggests an external origin f or the Hi 
gas (|Moran et al.ll2010l : see also IWang et al.ll2011^ . As can 
be seen, from a dynamical point of view GASS 35981 is not 
unusual, whereas GASS 3505 is still an outlier. 

We now investigate how the residuals of the BTF rela- 
tion depend on a few representative galaxy parameters. The 
left column of Figure[3]shows the BTF residuals for the sam- 
ple in Figure [2^ with respect to the best inverse fit obtained 
for the TF subset in Figure [2]3. In other words, residuals 
are computed as Log x — Log xat, where x is the measured 
rotational velocity and xat is the value expected from our 
best BTF relation for a galaxy with the same baryonic mass. 
Cyan filled circles indicate objects with inclination smaller 
than 40° . From top to bottom, residuals are plotted as func- 
tions of concentration index, inclination, NUV— r color, stel- 
lar mass surface density, gas fraction, and distance from the 
gas fraction plane mentioned above and described in Pa- 
per I (galaxies more Hl-rich than the average have positive 
distance), respectively. BTF residuals do not exhibit strong 
dependence on structural and star-forming galaxy proper- 
ties, except for a mild tendency toward increased scatter 
for more bulge-dominated, red galaxies, which largely disap- 
pears when the systems with low inclinations are removed 
from the sample. This is not surprising, as many attempts 
to identify a third parameter to minimize the scatter of the 



TF r e lation produced negative results (e.g., Courteau fc Rixl 
1 19991 : iPizagno et aLll2007i : iMever et al]|200a 



and references 

therein) . 

The results presented in this section show that the BTF 
is not very promising for our purpose of determining a re- 
lation between baryonic and dynamical mass that holds for 
all massive galaxies. Firstly, its application is restricted to 
galaxies with Hi detections and inclinations larger than 40° . 
Secondly, the bulge-dominated systems with Hi detections 
are not simply offset from the relation defined by the disk- 
dominated galaxies, but are also more scattered. As demon- 
strated in the following two sections, the BFJ relation does 
not suffer from these limitations, and can be generalized to 
include both disk-dominated and spheroidal systems. 

3.2 Baryonic Faber-Jackson Relation 

We now carry out a similar analysis for the BFJ relation. 
Since all the galaxies have a measurement of the stellar ve- 
locity dispersion from SDSS, we can use here the full sam- 
ple, but it is instructive to keep Hi detections and non- 
detections separated. In Figure |4] the baryonic mass is plot- 
ted as a function of the stellar velocity dispersion, and green 
upside-down triangles indicate galaxies that were not de- 
tected with Arecibo. The relation plotted for the full GASS 
sample (panel a) shows a clear segregation between Hi de- 
tections and non-detections, the former being offset towards 
lower values of velocity dispersion. The gas-rich elliptical 
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Figure 3. Residuals of baryonic TF (left column) and FJ (right column) relations plotted as functions of concentration index (row 
1), galaxy inclination (row 2), NUV-r color corrected for Galactic extinction only (row 3), stellar surface density (row 4), gas fraction 
(row 5), and distance from the gas fraction plane (see text; row 6). Filled cyan circles and green upside-down triangles indicate galaxies 
with inclinations smaller than 40° and Hi non-detections, respectively. For the top right panel, we show a linear fit to the data points 
(long-dashed line) and note its dispersion on the bottom right corner (see text). 
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Figure 4. Relation between baryonic mass and stellar velocity dispersion for our sample (a) and for selected subsets: (b) bulge-dominated 
systems with inclinations smaller than 70°, (c) galaxies with large inclinations, and (d) disk-dominated objects. Green upside-down 
triangles indicate galaxies that were not detected in Hi. The red star in (a) is GASS 3505; filled symbols in (c) indicate galaxies that 
are both bulge-dominated and with high inclinations, and blue filled circles in (d) are disk-dominated galaxies with large gas fractions 
{Mm/ Mir> 30%). In the top panels, a dotted (a) or dashed (b) line indicates the inverse fit to the data points (the scatter is noted in 
the bottom right corner). The dashed line in (b) is reproduced in panels (c) and (d) for reference. 



GASS 3505 (red star) is not unusual in terms of its stellar 
dispersion (GASS 35981 is not plotted because it has ct < 70 
km s"'^). We recover a tighter relation when we restrict our 
sample to the FJ subset (panel b), which includes bulge- 
dominated galaxies (i.e., objects with Rgo/R50> 2.8) with 
inclinations smaller than 70°. We excluded from the subset 
highly flattened systems because these are likely to host a 
signiflcant disk component, despite their high concentration 
index. The 26 bulge-dominated galaxies excluded by our in- 
clination cut are shown as flUed symbols in Figure These 
are mostly Hi detections, which already suggests that they 
are more likely to be disks than oblate spheroids. Their SDSS 
images (Figure [SJ conflrm that they are inclined disks, often 
with dust lanes running along their major axis. As demon- 
strated by the bottom panels of Figure 3] disk-dominated 
galaxies are systematically offset from the best fit relation 
obtained for the data points in panel (b), and they are 
mostly Hi detections. These objects have low stellar velocity 
dispersions for fixed baryonic mass, thus they increase the 



scatter of the FJ relation when they are plotted together 
with the bulge-dominated galaxies. 

Because FJ studies usually do not include the gas com- 
ponent, we check our results by considering the stellar FJ 
relation. In particular, we compared our FJ relation with 
the one obtained for a large sample of SDSS galaxies by 
iGallazzi et all (|2006l ). which is conveniently expressed in 
terms of stellar masses. The slope and scatter of our FJ re- 
lation (restricted to the FJ s ubset) are in re markable agree- 
ment with those published bv lGallazzi et al] (2006) (both re- 
lations have a scatter of 0.071 dex), w hereas the zero points 
differ by 0.17 dex in solar masses. The IGallazzi et al] (|2006l ) 
sample includes galaxies with concentration index > 2.8, 
redshift 0.005 < z < 0.22, and whose SDSS spectra have a 
median signal-to-noise per pixel greater than 20. Removing 
our inclination cut does not bring the zero points of the two 
relat ions into agreement (the offset is still 0.14 dex). How- 
ever, IGallazzi et al](|2005l ) note that their stellar masses are 
systematically larger than those derived bv lKauffmann et alj 



© 0000 RAS, MNRAS 000, 000-000 




The GALEX Areciho SDSS Survey. IV. 9 




Figure 5. SDSS postage stamps (1' size) for a subset of bulge-dominated galaxies with inclinations larger than 70° shown in Figure[4}; 
(filled symbols). The galaxies are labeled with their GASS identifier. Top three rows: galaxies with Hi detections, ordered by increasing 
inclination. Bottom row: Hi non-detections, also ordered by increasing inclination. 



o 




2.0 2.5 
Log (T [km s"'] 

Figure 6. Relation between bulge mass and velocity dispersion. 
Symbols and dotted line as in Figure |4^. 



l|2003ll (and adopted here) by ~0.1 dex, thus our results are 
completely consistent. Lastly, we notice that adding the gas 
to the stellar mass slightly increases the scatter of the rela- 
tion (from 0.071 dex for the FJ to 0.074 dex for the BFJ). 
This is because the most Hi-rich galaxies in our sample (see 
for instance the blue points in Fig. [4ji) become even fur- 



ther displaced from the Hi non-detections, which constitute 
the bulk of the FJ subset. A relation with similar scatter is 
obtained when, instead of restricting the sample to bulge- 
dominated systems, we consider only the bulge component 
itself. The mass of the bulge can be estimated from the cor- 
relation between SDSS concentration ind ex and bulge-to- 
total stellar mass ratio, B/T, presented bv lWeinmann et al.l 
l|2009l . see their Fig. 1), which is based on the 2D multicom- 
ponent decomposition analysis of iGadottH (|2009| ). Our bi- 
linear fit to their data yields B/T=(7?9o/ii50-1.920)/2.276. 
The relation between bulge mass and stellar velocity disper- 
sion plotted in Figure[S]has a scatter of 0.075 dex in velocity, 
comparable to that of t he FJ and BFJ relations app lied to 
the FJ subset (see also iGadotti &: Kauffmannll2009l ). Inci- 
dentally, it is somewhat surprising that we can recover such 
a tight relation given the uncertainties in our crude B/T 
estimate. 



We now focus again on the baryonic FJ relation pre- 
sented in Figure U and study its residuals. The BFJ residu- 
als are plotted as functions of Rgo/Rso, inclination, NUV— r 
color, stellar mass surface density, gas fraction and distance 
from gas fraction plane on the right column of Figure O 
As for the BTF, these are the residuals for the sample 
in Figure [4^ with respect to the best fit relation in Fig- 
ure 2)3, computed as described in § 13.11 (where x is now 
stellar velocity dispersion). As in all the figures in this pa- 
per, green upside-down triangles indicate Hi non-detections. 
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Interestingly, BFJ residuals show clear trends with concen- 
tration index, NUV— r color, /i*, and gas fraction, in the 
sense that more disky, gas-rich and star-forming galaxies 
have larger residuals. A similar trend is seen as a function 
of D„4000 index, which is an indicator of the age of the 
stellar population sampled by the SDSS 3"-diameter fiber 
(larger residuals are seen for smaller values of D„4000, not 
shown). It is well known that the FJ relation for ellipti- 
cal galaxie s is a projection of a mo r e general Fundamen - 
tal Plane l|Diorgovski fc David [l987l : iDressler et al.lll987l ). 
which is usually parametrized in terms of stellar velocity dis- 
persion, effective radius r^, and luminosity or surface bright- 
ness (central, or average within r^). Thus, one might expect 
to see a dependence of the FJ residuals on a third parame- 
ter related to r^. However, the trends in Figure |3] are not a 
consequence of the Fundamental Plane. Firstly, our sample 
does not include only elliptical galaxies and spheroids. As 
can be seen by inspecting Figure O the trends are driven by 
the disk-dominated, star-forming galaxies. Secondly, quan- 
tities such as Rgo/Rso, NUV— r color or gas fraction are not 
directly related to an effective radius. 

The comparison of the BTF and BFJ relations for 
our sample of massive galaxies illustrates several important 
points: (a) the BFJ is significantly less scattered than the 
BTF when the relations are applied to their ma:ximum sam- 
ple. When the two relations are applied to their respective 
"good", morphologically-pruned subsets, the scatter in both 
is almost identical, (b) the BFJ is insensitive to the inclina- 
tion problems that plague the BTF, which can be applied 
only to systems with inclination larger than 40°. Further- 
more, stellar dispersions are measured also for galaxies with- 
out Hi detections. Naturally, one could measure rotational 
velocities with other tracers and methods, e.g., with Ha ro- 
tation curves. However, these have their own sets of prob- 
lems (for instance. Ha emission is typically significantly less 
extended than Hi emission, thus it may not trace the full ro- 
tational velocity. Moreover, it is unclear at which spatial po- 
sition the rotational velocity should be measured. See, e.g., 
ICatinella et al.|[2007l ). and the inclination issue remains, (c) 
Most importantly, and contrary to the BTF case, the BFJ 
residuals show systematic trends with other galaxy proper- 
ties. As shown in Figures [3] and [J] disk-dominated galaxies 
do not form a scatter plot in the BFJ plane, but they sys- 
tematically deviate from the main relation defined by the 
bulge-dominated systems. 

One might wonder if, for disk-dominated, inclined 
galaxies, a is strongly affected by rotation. In other words, 
the BFJ for highly inclined disks (Figure HJ; and |3Jl) could 
just be a TF relation in dis guise. From the average rotation 
curves of disk galaxies of Catinella et al.l (|2006l ). we esti- 
mated that the typical rotational velocity reached at 1.5" 
by GASS objects is ~100 km s~^ (but could be up to twice 
as large for the most luminous and/or most distant galaxies 
in our sample j3- However, this does not take inclination and 



^ We converted SDSS model r magnit udes into Cousins / mag- 
nitudes following ICatinella et ahl | |200^) . and adopted Rgo,it the 
radius containing 90% of the Petrosian flux in the SDSS i-band, 
as a proxy for the optical radius -Ropt , which is the radius encom- 
passing 83% of the total /-band light of the galaxy. The average 
/-band absolute magnitude for our sample is —22.35 mag, and the 
average coverage of the SDSS fiber radius is 14% Ropt- 



12.0 



„ 1 1.5 




1 1.0 



+ 



10.5 



10.0 



9.5 

2.0 2.5 
Log (J^„„ [km s"'] 

Figure 7. Empirical baryonic FJ relation for massive galaxies, 
obtained by correcting the velocity dispersions for the systematic 
dependence on Rgg / R^q shown in Figure|3](top right panel). As in 
other figures, Hi non-detections are indicated as green triangles. 



bulge-to-disk ratio within the aperture into account. Most 
importantly, the fiber measurements are intensity-weighted, 
thus most of the disk stars contributing to the dispersion are 
those closer to the galaxy center, where Kot is negligible. 
Thus we conclude that a measurements are not significantly 
affected by contamination from disk stars in circular orbits. 

In summary, the SDSS stellar velocity dispersions pro- 
vide an estimate of dynamical mass that — at least in the 
stellar mass and redshift regime probed by GASS — is not 
only more generally applicable (i.e., not limited to inclined 
galaxies with Hi detections or with extended Ha rotation 
curves), but also less affected by measurement problems 
(which are largely responsible for the BTF outliers) com- 
pared to Vrot- Moreover, the fact that the disk-dominated 
galaxies follow a BFJ relation that is simply offset from 
that of the bulge-dominated systems suggests a simple way 
of bringing the two onto the same relation. This is accom- 
plished in the next section. 



3.3 A Generalized Baryonic Faber-Jackson 
Relation for All Massive Galaxies 

Spheroids and disk-dominated galaxies can be brought onto 
the same BFJ-like relation by correcting the velocity dis- 
persions for the trends observed in Figure [31 The tightest 
dependencies of the BFJ residuals are seen as a function of 
NUV— r color and concentration index. Although the corre- 
lation with NUV— r is slightly tighter (the dispersion of the 
linear fit is 0.073 dex), we decided to correct for the depen- 
dence on i?go/i?5o, which is a galaxy structural parameter. 
In fact, it is reasonable to expect larger corrections for more 
disk-dominated objects, where rotation is likely to be the 
main contributor to dynamical support. On the other hand, 
NUV— r color is linked to the star formation properties of 
the galaxy, thus the rationale for correcting velocity disper- 
sions based on this quantity is much less evident. Moreover, 
a correction based on NUV— r has the additional disadvan- 
tage of requiring UV photometry. 

The best linear fit to the data points in the top right 
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panel of Figure [3] is: 
Aa = 
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Figure 8. Comparison between data (circles and triangles, corre- 
sponding to Hi detections and non-detections, respectively) and 
simulations (red stars, see text): BFJ relation (a), BFJ residuals 
as function of R90/-R50 (b), and generalized BFJ relation (c). 



where the residuals Act = Log a — Log ctbfj are computed 
with respect to the best fit BFJ relation in Figure |4}3. We 
thus corrected the velocity dispersions Log a of all the galax- 
ies in our sample by subtracting the offset Act. The BFJ 
relation obtained using the corrected velocity dispersions is 
plotted in Figure [7] The inverse fit is almost indistinguish- 
able from that of Figure [Ifj, and has a scatter of 0.079 dex 
(the scatter in solar masses, obtained from a direct fit, is 
0.218 dex). This relation holds for all the galaxies in our 
sample, regardless of morphology, inclination or gas content, 
and is the main result of this work. 

We c ompared our resul t s with the simulations per- 
formed bv lScannapieco et al. I (|2009| . |2011| ). which follow the 
formation of eight Milky Way-mass halos in a A-cold dark 
matter cosmology, including baryonic physics (star forma- 
tion, metal cooling, chemical enrichment, multiphase gas, 
thermal feedback from supernovae). For these simulated 
galaxies we can measure stellar and cold gas masses, and 
stellar velocity dispersions at a given radius, and calculate 
concentration indices (from the ratio of the radii enclosing 
90% and 50% of the total luminosity; we use d r-band lu- 
mino sities computed from the dust-free Bruzu al fc Charlotl 
I2OO3I population synthesis models). We estimate total masses 
and luminosities within 10 kpc (the mean value of -R90 for 
the galaxies in the GASS sample), and obtain luminosity- 
weighted stellar velocity dispersions within 1.5 kpc (the 
physical size subtended by a 1.5" radius at z=0.05). Us- 
ing the concentration indices, we correct the stellar disper- 
sions according to Equation |4l The results are illustrated 
in Figure (HI where we compare the positions of the simu- 
lated galaxies (red stars) with those of the FJ subset on the 
BFJ plane (a). The residuals from the BFJ fit are plotted 
as a function of R90/R50 in (b), and the generalized BFJ 
relation is shown in (c). The simulated galaxies are all disk- 
dominated according to our definition (i.e., Rgo/Rso < 2.8) 
except one, and their stellar dispersions are systematically 
offset from the fit in (a) towards smaller values, as the disk 
galaxies in the GASS sample are. The correction applied to 
their stellar dispersions removes the offset, and the simu- 
lated galaxies lie on top of the relation in (c). Although it is 
based on eight halos only, the excellent agreement between 
simulations and data is encouraging. 

We have shown that disk-dominated galaxies are off- 
set in the BFJ plane, and that this is what allow us to ob- 
tain a tight baryonic mass- velocity relation that holds for all 
the galaxies in our sample. The reason why there is a disk 
BFJ relation for massive galaxies is that ct is proportional to 
Kot, and their ra tio is a function of galax y morphology. As 
demonstrated by ICourteau et al. I l|2007bl ). for a given Kot, 
earlier type galaxies have higher ct. The brightest, bulge- 
dominated galaxies (ellipticals, lenticulars and early-type 
spirals) lie on the Vrot = \/2ct relation expected for isother- 
mal stellar systems. Later-type spiral and dwarf galaxies, on 
the other hand, depart from the isothermal relation by an 
amount that depends on morphology or total light concen- 
tration (see their Fig. 1). We plot the relation between Vrot 
and CT for the Hl-detected GASS galaxies in Figure [9l Our 
sample shows two populations of outliers, characterized by 
too large or too low Kot/CT ratios compared to the rest of 
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Figure 9. Relation between rotational velocity and velocity dis- 
persion for GASS detections. Symbols are the same as those in 
Figure[2l Dashed and dotted lines indicate the 1:1 correlation and 
the Vrot = \/2o" relation, respectively. Error bars are plotted only 
for galaxies with inclination i smaller than 40° (cyan) or with 
Vrot < cr (red) that are not affected by beam confusion (magenta 
diamonds). Red symbols with cyan error bars are galaxies with 
both Vrot < CJ and i < 40° . 

the data set. These are the same outliers of the BTF rela- 
tion (see Fig. [2]); in particular, galaxies with large Kot/o" 
ratios all have small inclinations (< 40°, cyan), and objects 
with Vrot < c" (red) were mentioned in § 13.11 and are dis- 
cussed in the Appendix. To avoid crowding, we show error 
bars for these galaxies onljQ. Disregarding the outliers, one 
can see that the galaxies with larger stellar dispersions lie 
very close to the isothermal relation (indicated by a dotted 
line), whereas there is a tail of galaxies with higher Kot/o" 
ratios at lower dispersions. The galaxies in this tail are all 
disk dominated. Correcting the stellar dispersions according 
to Equation [4] effectively translates into accounting for this 
departure from isothermality for disk-dominated galaxies. 

3.4 Empirical Baryonic Mass - So. 5 Relation 

In the previous section, we obtained a relation between bary- 
onic mass and a measure of internal velocity that holds for 
our sample of massive galaxies, without any morphological 
pruning. We now ask whether we can find an even tighter 
relation by combining Vrot and a measurements. Both quan- 
tities measure the depth of the gravitational potential well 
of a galaxy, but they do so at different radii, namely those 
of the bulge and of the outer disk (as traced by the Hi gas, 
which typically extends well beyond the stellar disk; how- 
ever in very hi gh-density environments H i could be severely 
stripped; e.g., iGiovanelli fc Havned [l985l ) . In principle. Hi 
widths should provide the best measurement of the total 
mass (including dark) of a galaxy, regardless of the presence 

* The errors on rotational velocities are obtained from standard 
error propagation, taking into account the uncertainties in the 
observed velocity width and inclination (i.e., b/a axis ratio in 
Eq. 1) only. The errors on stellar velocity dispersions are from 
the SDSS (as noted in ii l2.2l the aperture correction is very small, 
thus we do not propagate its uncertainty). 



of a bulge component. However, as discussed above, reli- 
able Kot measurements are not possible for galaxies with 
small inclination to the line-of-sight — in the limit of a per- 
fectly face-on system, which has no measurable rotation, the 
velocity dispersion of the stellar component is a more use- 
ful quantity. Thus, we have two methods of estimating the 
dynamical mass of a galaxy that are affected by different 
systematics and limitations. Can we combine them in or- 
der to obtain a quantity that is, on average, more tightly 
correlated with the baryonic mass? An int eresting quantity 
in th is context is the So. 5 parameter (e.g.. ICovington et al.l 
|2010D : 

S0.5 = sj0.5V,l,+a^ (5) 

This expression has been used in several stud ies, with dif- 
ferent meanings for the dispersion component. iKassin et al.] 
studied the TF relation for an emission line-selected 
sample of 544 galaxies at 0.1 < z < 1.2, and demonstrated 
that a remarkably tighter relation is obtained when the So. 5 
parameter is adopted instead of the rotational velocity from 
the optical rotation curve. In their case, a is the disper- 
sion of the gas (but this may be contaminated by velocity 
gradients, see their section 4.2). Interestingly, their sample 
includes early to late spirals, irregular galaxies and merging 
systems, and their So.s-stellar mass relation for the lowest 
re dshift bin is in very g ood agreement with the FJ relation 
of iGallazzi et al.l l|2006l ) in term s of slope and z e ro po int. 

More relevant for our work, IZaritskv et al.l (|2008l ) used 
the So. 5 parameter (which they refer to as internal veloc- 
ity, V) with a measuring the stellar velocity dispersion of 
spheroidal galaxies. They showed that all classes of galax- 
ies lie on a two-dimensional surface, which is expressed as 
a linear combination of the logarithms of the effective ra- 
dius Tc, the internal velocity squared, the surface brightness 
within To, and the mass-to-light ratio within ro. Their sam- 
ple is a heterogeneous collection of 1925 spheroids and disk 
galaxies from existing data sets, which span the full range 
of galaxy types and luminosities. However, although they 
define V'^ = 0.5V,^ + cr^, they use either the circular velocity 
Vc for disk galaxies or the stellar velocity dispersion a for 
spheroids, but never add the two contributions. 

Another in t eresti ng analysis was carried out by 
ICovington et all (|2010l ). who studied the evolution of the 
stellar mass TF relation in a simulation of a disk merger. 
They argue that, since early- type galaxies are generally as- 
sumed to form through mergers of late- type systems, the 
scaling relations of early types should descend from those 
of late-type galaxies. Thus, they simulate the evolution of 
a galaxy merger, mimic observations of emission lines, and 
study how the kinematics of the system changes with time. 
The progenitors are two identical Sbc galaxies lying on the 
stellar mass TF relation, which merge and form a rotating 
elliptical galaxy. Intriguingly, they show that, while rota- 
tion is converted into stellar dispersion, the So. 5 parameter 
is approximately conserved — suggesting that So. 5 might be 
really tracing the mass distribution. 

The above studies motivated us to look at the relation 
between baryonic mass and So. 5 parameter for the galaxies 
in our sample, with an important caveat. The So. 5 param- 
eter defined by equation [5] has a physical meaning only if 
rotation and dispersion are associated to the same system. 
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Figure 10. So.5 relations for tiie full sample (a), for the subset of galaxies with inclination larger than 40° (b), and for the TF (c) and 
FJ (d) subsets. Hi non-detections are indicated as green triangles; the red star in (a) corresponds to GASS 3505. Dotted lines are inverse 
fits to the data points (the scatter is noted in the bottom right corner of each panel). 



For instance, for a rotating isothermal sphere the dynamical 
support is provided by both ordered and thermal motions, 
and the above combination is a direct result of the virial the- 
orem. For our galaxies, Viot measures the rotational velocity 
of the gas, and a is the velocity dispersion of the stars in the 
bulge, thus combining the two has no physical motivation. 
However, we are not trying to partition the dynamical sup- 
port of a galaxy between its bulge and disk components — 
as we already pointed out, Kot should account for the full 
dynamical support, regardless of the presence of a bulge. We 
wish instead to establish if we can average two measures of 
the potential well of a galaxy, which are based on different 
tracers, in order to obtain a quantity that, empirically, bet- 
ter correlates with the baryonic mass. After all, we expect 
the Hi width to be a more reliable measurement for inclined 
disk galaxies, and a to be better for elliptical or face-on, 
disk galaxies. Thus, the combination of the two might work 
better on average. 

The correlation between baryonic mass and S0.5 param- 
eter is plotted in Figure [TUl for the full sample (a), for the 
galaxies with inclinations larger than 40° (b), and for the 
TF and FJ subsets (c, d) defined in § [34] and § [SJ] The 



So.5 relation for the full sample has a scatter of 0.094 dex, 
very close to that of the BFJ relation for the full sample, but 
significantly larger than that of the generalized BFJ plotted 
in Figure [7| The S0.5 relation still suffers from the incli- 
nation problems that affect the BTF, although to a lesser 
degree. The strongest outliers seen in Figure llOb are re- 
moved by our 40° inclination cut (b,c). Notice that, when 
applied to disk-dominated galaxies with inclinations larger 
than 40°, the So. 5 relation has a scatter of only 0.064 dex, 
i.e., it is tighter than the BTF shown in Figure [2}3 (which 
has a scatter of 0.076 dex) — this is in fact the tightest rela- 
tion presented in this work (for comparison, the generalized 
BFJ relation restricted to the same sample has a scatter 
of 0.077 dex). Conversely, restricting the sample to bulge- 
dominated galaxies does not decrease the scatter of the S0.5 
relation. There is still a small offset between galaxies with 
and without Hi detections. Naturally, we do not have any 
information on the rotational velocity of the non-detections 
(some of which are disks close to edge-on view) — thus, for 
a fraction of these, So. 5 is certainly underestimated. Overall, 
the baryonic So. 5 relation is as good as our generalized BFJ, 
but only for galaxies with inclinations larger than 40°. 
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Figure 11. Relation between optical size and rotational velocity for all the galaxies with Hi detections (a) and for the TF subset (b). 
Symbols are the same as those in Figure [2] Dashed lines are bisector fits to the data points (the scatter is noted in the bottom right 
corner). 



4 SIZE- VELOCITY RELATIONS 

In this work we demonstrated that a simple correction ap- 
plied to the stellar velocity dispersions removes the offset 
between disk-dominated galaxies and spheroids seen in the 
BFJ relation, i.e., to first order, it removes its dependence on 
the internal structure of the galaxy. We further investigate 
this by considering the size-velocity relation. The existence 
of suc h a correlation is well e s tablished for both disk ga laxies 
fe.g.. ICourteau et al.|[2007al : lAvila-Reese et al.ll2008l : SSll 
and references therein) and ellipticals or bulges of early-type 
spirals (the Dn-cr relation, where Dn is the d iameter of the 
galaxy at a given surface br ightness level; e.g.. lDressler et al] 
Il987l : iBernardiet alll2002l ). Here, we compare this relation 
for the different velocity indicators and data subsets dis- 
cussed in this work. Specifically, we wish to test how the 
correction to the stellar velocity dispersion introduced in 
§ 13.31 affects the size- velocity relation. 

The sizes of disk galaxies are commonly estimated from 
exponential scale lengths, Vd- However, di sk scale lengths 
are n ot only problematic to measure (e.g. iGiovanelli et al.l 
Il994; SSll), but are also not meaningful for elliptical galax- 



Saintonge et al. I iipos) and SSll advocate the use of 



isophotal radii instead of disk scale lengths, as they yield 
scaling relations with significantly lower scatter. Isophotal 
radii are also used as size indicators for early-type galaxies 
in the Dn-cr relation. Thus, we adopt R25, which is half the 
25 mag arcsec"^ isophote diameter D25 measured by us on 
the SDSS g-band images, as a measure of size for all the 
galaxies in our sample. For comparison, R25 varies between 
2 and 7 rd, where rd is the SDSS exponential scale length 
in r-band, for the disk-dominated galaxies in our sample. 

The relation between size and rotational velocity (RV) 
is shown in Figure [TTK for all the galaxies with Hi detections. 
For this and the other relations presented in this section we 
performed bisector linear regressions instead of inverse fits, 
because it is less clear that the scatter is mostly confined to 
one variable. The scatter of R25 about the best fit is com- 
puted as before, i.e., by applying Tukey's bi- weight, and is 
noted in the bottom right corner of the plot (in dex of kpc) . 



As in Figures [2] and (5] cyan and red circles indicate galax- 
ies with inclinations smaller than 40° and with Vrot/<y < 1, 
respectively. Some of the main outliers of this relation are 
also outliers for the BTF, but there is not a 1-to-l corre- 
spondence. For instance, the very Hl-rich spiral GASS 35981 
(marked as a blue star) is an outlier on this plot, but not 
on the BTF, whereas the gas-rich elliptical GASS 3505 (red 
star) is an outlier for both relations. 

The RV relation is significantly weaker and more scat- 
tered than the BTF, as known from other studies t hat use 
disk s cale lengths as size indicators (e.g., Courtoa u et al.l 
l2007al report a Pearson correlation coefficient of r ~ 0.65 
and a scatter of 0.33 in log R). However, as mentioned above, 
SSll obtain a very tight correlation using isophotal radii in- 
stead of disk scale lengths (r = 0.84 and a scatter of 0.11 in 
log R; in the notation adopted by these authors, the scatter 
of Courteau's relation becomes 0.165 dex). When restricted 
to the TF subset (Fig. Illb ). the scatter of our relation 
(0.15 d ex) is intermediate between those of lCourteau et al] 
(2007a) and SSll samples. Although we use a similar size 
indicator as SSll, their sample includes only late-type spi- 
rals, thus the larger scatter of our relation is most likely due 
to the broader morphological mix of the galaxies in our data 
set. As for the BTF and BFJ relations, we plot the resid- 
uals of the RV relation as a function of several quantities 
in Figure [12] (left panels) . There is a dependence of the RV 
residuals on inclination and, to a lesser extent, on /x*. 

Figure [T3] shows the relation between R25 and stellar 
velocity dispersion a (Ra) for the full sample (a) and for 
the FJ subset (b). As for the BFJ relation, galaxies with Hi 
detections are displaced from the non-detections (green tri- 
angles), and a clear trend is observed when the Ro" residuals 
are plotted as a function of concentration index (Fig. 1121 
top right panel). Indeed, Figure \T2\ shows that the Ra resid- 
uals behave similarly to the BFJ ones (see Fig. [Sjl: disk- 
dominated galaxies systematically deviate from the best fit 
relation obtained for the FJ subset, which is dominated by 
spheroids. Aside from systematic trends, which are stronger 
for the R(j relation, the scatter of the RV and Rct relations 



© 0000 RAS, MNRAS 000, 000-000 



The GALEX Areciho SDSS Survey. IV. 15 




1.5 2.0 2.5 3.0 3.5 4.0 1.5 2.0 2.5 3.0 3.5 4.0 




20 40 60 80 20 40 60 80 

inci [deg] inci [deg] 



1.0 




0.5 




0.0 




0.5 




1.0 





2 3 4 5 

NUV-r [mag] 



b 

en 



1.0 

0.5 

0.0 

-0.5 
-1.0 



3 4 5 

NUV-r [mog] 




8.5 9.0 9.5 

Log jj., [Mq kpc"^] 



10.0 




8.5 9.0 
Log /i, [Mq kpc"^] 



10.0 




-2.0 -1.5 -1.0 -0.5 0.0 -2.0 -1.5 -1.0 -0.5 0.0 

Log Mhi/M. Log M„/M, 




-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 



Figure 12. Residuals of RV (left column) and Rct (right column) relations plotted as functions of concentration index (row 1), galaxy 
inclination (row 2), NUV— r color (row 3), stellar mass surface density (row 4), gas fraction (row 5), and distance from the gas fraction 
plane (row 6). The residuals are computed from the orthogonal distances to the bisector fits shown in Figurcs lllb and ll3b . respectively. 
Cyan and green symbols indicate galaxies with inclinations smaller than 40° and Hi non-detections, respectively. 
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Figure 14. Same as Figure [T3b for the corrected velocity disper- 
sion. 

is very similar, for both maximum samples and pruned sub- 
sets. 

Lastly, correcting the velocity dispersions according to 
Equation [4] yields a slightly tighter size-velocity relation 
(Fig. I14|) . Most importantly, our correction removes the off- 
set between Hi detections and non-detections present in Fig- 
ure 113b , (the remaining few Hi detections lying above the 
relation are highly inclined galaxies, for which R25 , which is 
not corrected for inclination, is likely to be overestimated). 



5 DISCUSSION AND CONCLUSIONS 

The main result of this work is the existence of a tight re- 
lation between baryonic mass and velocity that is indepen- 
dent of galaxy morphology. To first order, we can remove 
the dependence of dynamical scaling relations on the inter- 
nal structure of the galaxy by applying a simple correction 
to the measured stellar velocity dispersions, which depends 
only on the concentration index Rgo/Rso (Eq.[4l). The corre- 



lation between baryonic mass and corrected a thus obtained 
has a scatter of only 0.08 dex (Fig. [7] and Table [T]). It is en- 
couraging that our correction removes the offset between 
disks and spheroids also in the size-velocity relation (ex- 
pressed in terms of R25 versus stellar dispersion). We tested 
if an even tighter baryonic mass-velocity relation could be 
obtained by using the ^o.s parameter, which combines rota- 
tional velocity and stellar dispersion. We found no improve- 
ment, because the So. 5 relation still suffers from the problem 
of low inclination galaxies that affects the BTF, to a lesser 
extent. Despite being measured at a smaller spatial scale 
than the Hi width, the stellar velocity dispersion turns out 
to be a better tracer of mass, at least for the massive galax- 
ies in our sample. This is surprising, as one would expect the 
rotational velocity to provide a more reliable measurement 
of mass for galaxies with Hi detections, the vast majority 
of which are rotation dominated (i.e., 216/228 galaxies have 
Kot/o" > 1), regardless of the pre sence of a bulge. An d yet, 
the same result was obtained by Neistein et al.l (| 19991 ) for a 
sample of SO galaxies with inclinations ~35°— 60°, which are 
also, overall, rotation-dominated: the central stellar velocity 
dispersion is a better predictor of /-band luminosity than the 
circular speed at 2-3 exponential disk scale lengths, which 
they carefully measured from long-slit optical absorption- 
line spectra. 

As we already pointed out, the main limitation with 
rotational velocities is observational, not intrinsic. We are 
simply not able to reliably measure circular speeds from 
line-of-sight velocities, except for well-selected samples of 
undisturbed late-type, inclined spirals (i.e., the typical TF 
samples), for which the deprojection to edge-on view does 
not introduce very large uncertainties. As for the velocity 
dispersions, which are measured through 3"-diameter fibers, 
there might be a concern about contamination from disk 
stars in circular orbits, especially for the most edge-on, disk- 
dominated galaxies. However, we have argued that this ef- 
fect is negligible (§ 13.21) . T his agr ees with the conclusion 
reached bv lCourteau et al.l l|2007bl ). who reported that the 
contamination is small (< 5%, based on simulations) for 
Milky Way-type galaxies with pressure-supported bulges. 

We argued that the reason why disk-dominated galax- 
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ies are simply offset from the spheroids on the BFJ plane 
(which is the key point that allows us to bring all massive 
galaxies onto the same baryonic mass-velocity relation) is 
that a is proportional to Vrot, and their ratio is a function 
of galaxy morphology. Although initial work indicated that 
spirals and ellipticals follo w the same tight Vrot-cr correlation 
for g > 80 km s~^ (e.g.. lFerrare5 12002| : [Baes et al.ll2003l : 
iPizzella et al. more recent analyses demonstrate that 

the relation between Kot and a is not un iversal, but de- 
pends o n morphology. As mentioned in i) l3.3|[Courteau et al.l 
l|2007bl ) show that the brightest, bulge-dominated galaxies 
lie on the K-ot = \/2(t relation expected for isothermal stellar 
systems, whereas later-type spirals and dwarfs are offset by 
an amount that depends on morphology or total light con- 
centration. We have shown that correcting the stellar dis- 
persions according to Equation |4] effectively translates into 
accounting for this departure from isothermality for disk- 
dominated galaxies. ICourteau stall (|2007bl l noted that, de- 
spite the fact that a detailed understanding of what sets 
the relation between Kot, cr and concentration index is still 
missing, one could use this relation to empirically reduce the 
scatter of scaling relations that involve dynamical parame- 
ters, such as the TF or FJ. We showed that the existence of 
such relation allows us to do more than that — we obtained 
a generalized BFJ relation that holds for all the massive 
galaxies in our sample, with a scatter (0.079 dex) that is 
as small as that of the BTF and BFJ relations applied to 
their pruned su bsets (i.e. 0.076 and 0.074 dex , respectively). 
For co mparison, I Avila- Reese et al] l|2008h and lGallazzi et all 
lj2006t ) report a scatter of 0.06 dex for the BTF and 0.071 
dex for the stellar FJ relations, respectively. 

The implications of our generalized BFJ relation for ex- 
tragalactic studies are very promising. Because it holds for 
all massive galaxies in our sample regardless of morphol- 
ogy, this relation appears to provide a more fundamental 
link between dark matter halo mass and baryonic content 
than that obtained using rotational velocities or velocity dis- 
persions. As such, it gives more fundamental constraints to 
galaxy formation models than the TF or FJ/FP relations. 
Also this is the reference baryonic mass-internal velocity re- 
lation that higher redshift studies, which naturally target 
the most massive galaxies, should compare with. This rela- 
tion is more resilient to systematic effects than the TF, and, 
contrary to both TF and FJ/FP relations, does not require 
any morphological pruning — a significant advantage since 
accurate morphological classifications are difficult to obtain 
for large samples beyond the very local Universe. 

As pointed out throughout this work, our results are 
based on a representative sample of galaxies with stellar 
masses larger than 10^° M0. It remains to be established 
how far down in stellar mass these results can be extrapo- 
lated. It would be beneficial to investigate whether a BFJ 
relation for disk galaxies holds down to low baryonic masses 
(where the gas contribution is more important), and with 
similarly low scatter, based on a representative and homo- 
geneous sample such as GASS. 
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APPENDIX: GALAXIES WITH Vrot < a 

We discuss here in more detail the 12 galaxies with Hi ro- 
tational velocities smaller than their stellar velocity disper- 
sions, which are outliers of both BTF and Kot versus a 
relations. These objects are indicated by red symbols in Fig- 
ures [2] and |9l While the high Viot /<y outliers are all galax- 
ies with small inclinations, for which the rotational veloci- 
ties are likely to be overestimated, the low Vrot/o outliers 
are po tentially more interesting. As mentioned in § 13.11 IHoI 
|2003) identifies a population of galaxies with unusually 
small Viot/iy ratios, which are outliers in his TF relation, 
and argues that these systems must have experienced gas 
accretion. Figure [151 shows SDSS images and Hl-line profiles 
for 9 outliers detected by GASS (the other galaxies were ei- 
ther observed by ALFALFA or included in the SOS archive). 
Two of these galaxies, GASS 30175 and 31156, are outliers 
because their SDSS inclinations are incorrect. They both 
have strong bars which dominate the SDSS fits, thus overes- 
timating the inclination (54° and 61° for GASS 30175 and 
31156; we measured 37° and 27° respectively). The same is 
true for GASS 24236 (AGC 220363, detected by ALFALFA 
and not shown here), a nearly face-on barred galaxy with 
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Figure 15. SDSS postage stamps (1' size) and Hi spectra of subset of outliers marked in red in Figures [2] and |9] The GASS identifier 
of each galaxy is indicated on the top right corner of its spectrum. A dotted line and two dashes indicate the heliocentric velocity 
corresponding to the SDSS redshift and the two peaks used for width measurement, respectively. First row and GASS 31156: GASS DRl 
(Paper 1); other objects: GASS DR2 (Catinella et al., in preparation). 



SDSS inclination of 46° (we measured 17°). Of the remain- 
ing outliers, one is clearly morphologically disturbed (GASS 
57099, bottom row). Another galaxy, GAS S 4030, was ob- 
serve d as part of the COLDGASS survey l|Saintonge et al.l 
I2OIII . see their Figure Al, row 6). Its (uncorrected) veloc- 
ity width measured from the CO (1-0) spectrum is 137 ±11 
km s"'^, which is exactly the same as the (uncorrected) Hi 
width. Since the CO (1-0) emission typically originates from 
the central regions of a galaxy, perhaps the fact that CO 
and Hi widths are identical means that the Hi gas in this 
galaxy is not extended enough to sample the flat part of 
the rotation curve. This is not particularly unusual, as there 
are other examples of early-type spirals with truncated Hi 



disks (e.g. NGC 3623, where the Hi emission does not ex- 
tend beyond the stellar disk; Hogg et al. 2001) or unexpect- 
edly narrow Hi profiles (e.g. NGC 5854; Hayncs ct al. 200(i). 
Overall, even acknowledging that some of these profiles are 
not very high signal-to-noise detections, it is striking that 
most are asymmetric, suggesting disturbances in the distri- 
bution and/or kinematics of the Hi gas. We will investigate 
asymmetries of Hi profiles in more detail in a future work. 
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